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Abstract 


This study analyzes the relationship between land use changes over different time periods and land cover metrics using 
data from the Barandozchay watershed in West Azerbaijan Province, Iran. To achieve this goal, Sentinel-2 satellite 
images for the years 2016 and 2022 were obtained from the European Union's Copernicus website. Pre-processing 
techniques were applied using various software tools, and the images were classified through knowledge-based and 
object-oriented methods in the eCognition software environment. Land use maps were generated, and land cover 
metrics were quantified and calculated at landscape and class levels using Fragstats 8.2 software. 

The results at the class level showed that the split index (SPLIT) had the highest value for the blue area, while the 
cohesion index (COHESION) had the lowest value in both years. The grassland class exhibited the highest SHAPE- 
MN and LSI indices, indicating patch disorder. If protective measures are not implemented, further degradation of 
grasslands in the Barandozchay watershed could occur in the coming years. 

At the landscape level, the results revealed a decreasing trend in most indices related to watershed connectivity. The 
ENN MN index, representing the mean Euclidean distance between nearest neighbors, the CONTAG index, 
representing contagion, and the COHESION index, representing patch cohesion, all showed a decline from 2016 to 
2022. Conversely, the LSI and NP indices displayed an increasing trend, indicating greater irregularity within the area. 
The findings also showed that the blue area experienced the highest level of fragmentation, whereas dry farming 
exhibited the least. Additionally, grassland patches demonstrated the highest connectivity, while dry farming patches 
had the lowest. At the landscape level, the results indicated a uniform distribution of patches in the Barandozchay 
watershed. The increasing number of patches, extended margins, and reduced patch size led to fragmentation and the 
formation of smaller, isolated patches.These outcomes highlight the ongoing land use changes and their impact on the 
structure and connectivity of the landscape. They provide valuable insights for land use management and conservation 
efforts in the study area, ensuring the sustainability of natural resources for future generations. 
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حکیده 


هدف از پژوهش حاضر. تحلیل ارتباط بین تغییرات کاربری اراضی در دوره‌های زمانی مختلف با استفاده از سنجه‌های 
سیمای سرزمین در حوضه آبریز باراندوزچای در آذربایجان غربی است. بدین منظور. تصاویر ماهواره‌ای سنتینل-۲ 
(Sentinel2)‏ مربوط به سال‌های ۱۳۹۵ و ۱۶۰۱ از Cole‏ کبورنیک اتحادیه ارویا دریافت گردید. سپس روش‌های 
پیش‌پردازش در lume‏ نرم‌افزارهای مختلف اعمال و تصاویر مربوطه به de‏ نرم‌افزار 0 ارسال شد. برای 
سنجه‌های سیمای سرزمین از نرم‌افزار 2 Fragstats‏ در دو سطح سیما YY)‏ سنجه) و کلاس (۱۳ سنجه) کمی‌سازی 
و محاسبه شدند. در سطح کلاس پهنه آبی. نتایج نشان می‌دهند در هر دو سال. بیشترین مقدار سنجه تکه‌شد گی 
(SPLIT)‏ و کمترین مقدار سنجه یکپارچگی (COHESION)‏ به دست آمده است. کلاس مرتع دارای بیش‌ترین 
مقدار در شاخحص (SHAPE-MN)‏ و (LSD‏ داشته است که این نشان‌دهنده بی‌نظمی در لکه‌های مرتع است. 
شاخص‌های (CONTAG) (ENN_MN)‏ و (COHESION)‏ که به ترتیب نشان‌دهنده‌ی میانگین فاصله 
نزدیک‌ترین همسایه اقلیدسی» شاخحص سرایت و شاخص پیوستگی لکه روند کاهشی را از سال ۱۳۹۵ نسبت به سال 
٤١١‏ تجربه کرده Eal‏ شاخص‌های شاخحص شکل لکه (LSI)‏ تعداد لکه ) بی نظمی را نشال می‌دهد که روند 
افزایشی داشته است و نشان‌دهنده افزایش بی‌نظمی در منطقه موردمطالعه است. نتایج حاصل از سطح سیما نشان داد 
که پراکنش لکه‌ها در حوزه آبخیز باراندوزچای CHES‏ است و افزایش تعداد لکه‌هاء افزايش حاشیه و کوچک شدن 
لکه‌ها باعث ES‏ و درنتیجه ایجاد لکه‌های جدا از هم در حوزه x‏ باراندوزچای می‌شود. این نتایج می‌تواند 
در مدیریت کاربری اراضی و حفاظت منطقه موردمطالعه برای نسل‌های آتی مفید باشد. 


Email: j.behmanesh ۵ urmia.ac.ir ۰۹۱۶۱۶۰۲۷۲ نویسنده مسئول:‎ . ` 


کوهانی و بهمنش. تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین e‏ ۳۶۹ 


کلیدواژه‌ها: تغییرکاربری اراضی. Ti tk‏ چشم‌انداز 3 تصاویر ماهواره‌ای سنتینل -۰۲ طبقه‌بندی شیء گرا 3 تحلیل 


۱- مقدمه 

تغییر کاربری زمین یکی از مهم‌ترین عوامل برای درک تغییر محیطی در تمام مقیاس‌های زمانی-مکانی 
است )2015 (Kumar, Bhaskar & Padmakumari,‏ تغییر کاربری زمین با عامل‌هایی مانند افزایش ارتفاع 
در حوضه‌های آبریز: افزایش ارتفاع آب بر حجم رواناب در رودخانه و افزايش ظرفیت سیلابی 
رودخانه‌ها تأثیر خواهد داشت(2017 (Idham et al.,‏ . یکی از مهم‌ترین پارامترهای تصمیم گیری برای 
برنام‌ریزی زمین در مقیاس جهانی و ناحیه‌ای آگاهی از تغییر کاربری زمین است. پیشرفت در درک 
حالت‌های ممکن تغییر کاربری زمین می‌تواند در شناسایی و تشخیص فرآیندهای زمین (سیلاب و ...) در 
آینده مناسب باشتد(2019 .(Armenteras, Murcia, González, Barón & Arias,‏ 

مهم است که چنین تغییراتی به‌طور دقیق ثبت شود تا فرآیندهای انسانی در کار بهتر درک شود. کاربری 
اراضی و پوشش اراضی ازجمله فرآیندهایی طبیعی و انسانی است که تغییرات در OF‏ خیلی سریع است 
.(Dezhbani, Mustafazadeh, Esmaaliouri & Nazila, 2023)‏ فناوری سنجش از oly‏ دور راهکاری ارزشمند 
در تهیه نقشه تغییرات کاربری اراضی است که امروزه در جهت شناسایی. پایش. تشخیص و پهنه‌بندی منابع 
طبیعی» به‌ویژه در مورد روند تهیه نقشه کاربری اراضی» استفاده می‌شود. این روش به دلیل هزینه‌بر و زمان‌بر 
بودن عملیات زمینی؛ جایگزین مناسبی برای تهیه اطلاعات دقیق و در سطح وسیع می‌باشد ( Nazar Neghad,‏ 
(Hosseine & Mostafazadeh, 2020‏ برای استخراج اطلاعات کاربردی طبقه‌بندی تصاویر ماهواره‌ای یکی از 
VE‏ روش‌ها محسوب می‌شود )2010 (Blaschke,‏ 

در حال حاضر طبقه‌بندی poled‏ رقومی ماهواره‌ای با سه روش کلی پردازش پیکسل‌پایه. شیءگرا و 
آموزش گرای انجام می‌شود )2018 f „=l .) Ghafari, Moradi & Modares,‏ روش جدید شیءگرا برای 
تجزیه‌وتحلیل تصویر استفاده می‌شود که قطعه‌بندی سطح پایین و بدون دانش را با روش‌های طبقه‌بندی فازی 
سطح بالا و مبتنی بر دانش ادغام می‌کند. این روش جدید از طریق نرم‌افزار eCognition‏ که شامل محیطی 


شیءگرا برای طبقه‌بندی تصاویر ماهواره‌ای op‏ طراحی گردید(2008 (Tzotsos & Argialas,‏ . در روش 


i ۳۷۰‏ 4 جغرافیا و مخاطرات محیطی. جلد ۰۱۳ شماره ۳ پائیز ۱۶۰۳ 


طبقه‌بندی شیءگرا اطلاعات eib‏ شکل و بافت‌ها را به‌طور هم‌زمان در نظر می‌شوند. در این روش 
خطاهای متعددی نسبت به روش‌های پیکسلی قابل‌انتظار است. در روش شیءگرا مجموعه‌ای پیکسل‌های 
مشابه تحت یک عنوان شناسایی و سپس عملیات طبقه‌بندی و تقسیم‌بندی اعمال می‌شوند ( Nazar Neghad‏ 
(et al, 2020‏ در تجزیه‌وتحلیل تغییرات کاربری qu‏ از عوامل طبیعی و اجتماعی_اقتصادی نیز استفاده 
می‌شود(2006 cal ab (Aunap, Uuemaa, Roosaare & Mander,‏ کاربری‌های مختلف زمین بر خدمات 
بوم‌سازگان‌ها (ES)‏ از طریق سه روش اصلی شامل الگوهای کاربری زمین. شیوه‌های استفاده از زمین و 
شدت‌های مختلف استفاده از زمین نمود پیدا می‌کنند ( Uuemaa, Antrop, Roosaare, Marja & Mander,‏ 
9 مطالعه تغییرات کاربری اراضی زمینه hel‏ پژوهش‌های سیمای سرزمین است که یکی از مسائل مهم 
در تغییرات بوم‌شناختی جهانی است(2023 (Dezhbani et al.,‏ نرم‌افزار Fragstats‏ یکی از ابزارهای محاسبه 
سنجه‌های سیمای سرزمین؛ می‌باشد(2002 (MeGarigal,‏ برای ارزیابی روند تغیبرات کاربری اراضی می‌توان 
از سنجه‌های سیمای سرزمین که یک رویکرد تلفیقی با سنجش ازدور است. استفاده کرد. این سنجه‌ها درواقع 
شاخص‌هایی هستند که Ga Shy‏ اجزای ساختار فضایی سیمای سرزمین مثل هندسه شکل. ماهیت 
پراکنش و توزیع این اجزا را به‌صورت کمی قابل‌مقایسه می‌سازند(  Kakehmami, Moamari, Qurbani‏ 
.(Ghaffari, 0‏ سیمای سرزمین شامل محدوده‌ای است که Lebe‏ زیست‌شناختی و فیزیکی. لکه‌ها و 
زیستگاه‌ها را در بر می‌گیرد. سنجه‌های سیمای سرزمین در سه سطح کلاس. سیما و aS‏ محاسبه می‌شوند. 
سنجه‌های سطح کلاس برای یک نوع کاربری» سنجه‌های سطح لکه ویژگی مکانی لکه‌ها و سنجه‌های سیمای 
سرزمین برآیند ویژگی‌های لکه‌ها و طبقه‌ها هستند(2018 (Nasiri & Darvishsefat,‏ . بسیاری از معیارهای 
الگوی سیمای سرزمین نسبت به تغییرات فضایی وضوح داده‌ها یا مساحت (وسعت) سیمای سرزمین حساس 
هستند و همبستگی‌های متعددی بین شاخحص‌های سیمای سرزمین رخ می‌دهد )2009 .(Uuemaa et al.,‏ 

لام و همکاران )2018 (Lam, Cheng, Zou & Cai,‏ به بررسی اثر تکه‌شدگی (SPLIT)‏ سیمای سرزمین بر از 
دست دادن زمین در ساحل Louisiana‏ در ایالت جنوبی آمریکا پرداختند. در این پژوهش از تصاویر لندست 


TM+‏ سال‌های 3441 تا ۲۰۱۰ و از ابعاد فراکتال منطقه تگی فضایی ' استفاده شد. نتایج OUS‏ داد که 
ی وارد و w^‏ یی e‏ 


' Spatial autocorrelation 
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تکه‌شدگی سیمای سرزمین منجر به از دست دادن بیش‌تر زمین می‌شود. Up‏ و ایسوفی (Berila & Isufi,‏ 
)2021 در مطالعه‌ای در منطقه کوزوو پراکندگی شهرها با استفاده از GIS‏ سنجش‌ازدور و سنجه‌های سیمای 
سرزمین را بررسی نمودند. نتایج نشان داد در طول دوره ۲۰۲۰-۲۰۰۰ در شهر پریشتینا ' مساحت مناطق 
شهری به میزان ۱۳/۶۷ کیلومترمربع افزایش abl,‏ است. در ایران نیز غفاری و همکاران ) Ghafari et al.,‏ 
2018( در دشت‌های اصفهان- برخوان نجف‌آباد و چادگان با استفاده از تصاویر لندست TM‏ سال ۱۹۸۵ و 
JLOLI‏ ۲۰۱۵ تغییرات کاربری اراضی را با استفاده از روش طبقه‌بندی شیءگرا و پیکسل ah‏ مقایسه کرد. 
نتایج نشان داد طبقه‌بندی شیءگرا با دقت بالای ٩۰‏ درصد نسبت به طبقه‌بندی پیکسل‌پایه. دقیق‌تر بوده است. 
نظرنژاد و همکاران(2020 (Nazar Neghad et al.,‏ به ارزیابی تغییرات پیوستگی و الگوی کاربری اراضی با 
استفاده از سنجه‌های سیمای سرزمین در حوزه آبخیز زولاچای. سلماس پرداختند. نتایج نشان‌دهنده افزایش 
سنجه تعداد لکه (NP)‏ و افزایش تراکم لکه (PD)‏ بود که این امر افزایش ناهمگنی در کاربری اراضی را نشان 
می‌دهد. همچنین اسفندیاری و همکاران(2021 (Esfandiari darabad, Hamzaei, Alaei & Mostafazadeh,‏ در 
ارزیابی تغییرات مکانی سنجه‌های سیمای سرزمین در پوشش گیاهی جوامع حاشیه متأثر از نوع کاربری 
اراضی در رودخانه قره‌سو استان اردبیل نتیجه گرفتند که تفاوت معنی‌داری بین کاربری باغ با سایر کاربری‌ها 
در سطح لکه و سیما وجود دارد. در مطالعه‌ی کیانی و فقهی )2015 (Kiyani & Feghhi,‏ بررسی ساختار 
کاربری اراضی حوزه آبخیز سفیدرود با استفاده از سنجه‌های بوم‌شناسی سیمای سرزمین بررسی شد. در این 
پژوهش, از نرم‌افزار ArcGIS‏ برای طبقه‌بندی مجدد نقشه کاربری اراضی حوزه آبخیز سفیدرود استفاده شد. 
سپس لایه‌ها به فرمت Grid‏ تبدیل و با استفاده از نرم‌افزار Fragstats‏ « سنجه‌های مربوطه استخراج شدند. نتایج 
نشان داد که بیشترین گسستگی مربوط به کاربری کشاورزی بوده و کمترین گسستگی مربوط به کاربری آب 
است. در مطالعه‌ی دژبانی و همکاران(2023 (Dezhbani et al.,‏ ارتباط بین تغییرات کاربری اراضی در آبریز 
کوزه‌تبراقی با استفاده از سنجه‌های سیمای سرزمین در دوره‌های زمانی مختلف بررسی شد. این مطالعه از 
تصاویر ماهواره‌ای لندست برای سال‌های ۱۳۷۹ ۱۳۸۹ و ۱۶۰۰ استفاده کرد. سنجه‌های سیمای سرزمین در 


دو سطح سیما و کلاس محاسبه شدند. نتایج نشان داد که در سطح کلاس» زراعت ديم و زراعت ual‏ در 


! Prishtina 


۳۷۲ نشریه جغرافیا و مخاطرات محیطی, جلد OY‏ شماره Y‏ پائیز ۱۶۰۳ 


سنجه‌های مختلف بیشترین مقدار را داشته‌اند. در سطح lage‏ شاخحص تکه‌شدگی در سال ۱۳۸۹ کاهش 
جزئی داشته و سپس در سال ۱۶۰۰ افزايش پیدا کرده است. شاخحص بزرگ‌ترین لکه نیز در سال ۱۶۰۰ 
کاهش قابل‌توجهی داشته است.بر اساس پژوهش‌های انجام IEE edd‏ بر لزوم انجام مطالعات گسترده‌تر و 
جامع تر در زمینه ارزیابی تغیبرات کاربری اراضی و نیز تحلیل ارتباط بین آن‌ها و تغییرات سنجه‌های سیمای 
سرزمین دارد. بر همین اساس» در پژوهش حاضر بررسی تغییرات کاربری اراضی در دوره‌های زمانی ۱۳۹۵ 
و ۱۶۰۱ با استفاده از سنجه‌های سیمای سرزمین برای حوزه آبخیز باراندوزچای در استان آذربایجان غربی 
انجام شد تا با اطلاع داشتن از تغییرات هر یک از کاربری‌هاء اقدامات و تصمیمات مدیریتی مناسبی برای 


استفاده درست از صورت گیرد. 


7- مواد و روش‌ها 


۱-۴- منطقه موردمطالعه 


حوزه آبخیز باراندوز چای Bly‏ در استان آذربایجان غربی است. گستردگی جغرافیایی این حوضه از ۰5۳۳۷ 
تا "۳۷۳۲۹ عرض‌های شمالی و "10*46 تا "10*۱۶ طول‌های شرقی واقع گردیده است. رودخانه باراندوزچای 
از دو رشته ارتفاعات مرزی olal‏ و ترکیه و ایران و عراق به نام‌های جمال‌الدین و ارتفاعات ککوداغ و الوق 
سرچشمه می گیرد و در بستر عمیقی با شیب تند > OS‏ کرده پس از دریافت شاخه‌هایی چند و طی مسافتی 
وارد جلگه شده دشت بیل را Solel‏ می کند و سپس از Goll‏ زیوه گذشته به روستای هفتوان باراندوز می‌رسد 
و از این ناحیه به بعد به نام باراندوزچای نامیده می‌شود. رودخانه باراندوز پس از طی مسافتی و آبیاری کردن 
زمین‌های متعدد روستایی. در چندین شاخه که شعبه hel‏ آن بابارود است به دریاچه اورمیه می‌ریزد. این 
حوضه در محدوده شهرستان‌های ارومیه و نقده قرار دارد. مساحت حوضه gSde‏ ۱۳۱۱/۵۶ کیلومترمربع 
است. آب‌های حاصل از این حوضه با جهت غربی - شرقی, درنهایت وارد دریاچه ارومیه می‌شوند. حوضه 
مذکور اغلب کوهستانی بوده و تنها بخش‌های شمال شرقی و مرکزی دارای ارتفاع پایین بوده و عمدتاً از 
دشت‌های حاصلخیز تشکیل شده است. این امر سبب شده تا در این قسمت‌هاء. پوشش گیاهی وسیعی باغات 


کشاورزی وجود داشته باشد. این حوضه عمدتاً از سنگ‌های «ubl‏ دولومیت‌ها. شیل‌ها. کنگلومراها و در 


کوهانی و بهمنش, تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... YWO‏ 


دشت‌ها از رسوبات کواترنری تشکیل شده است. شرایط ناهمواری» زمین‌شناسی. ژئومورفولوژی, فعالیت‌های 
انسانی و پوشش گیاهی شرایط مناسبی را برای وقوع پدیده‌های ژئومورفولوژیکی مختلف ایجاد کرده است. 
از لحاظ ویژگی‌های فیزیوگرافیکی حوضه باراندوزچای, دارای ۶ واحد GIS‏ بابارود. قاسملو» دیزج و 
هاشم‌آباد می‌باشد. بابارود با مساحت ۱۱2۹/۱۵ کیلومترمربع. بیشترین وسعت و قاسملو با مساحت ۳۳۱/۷۶ 
کیلومترمربع کمترین وسعت را به خود اختصاص داده است. بیشترین شیب متوسط وزنی با مقدار ۳۲/۹۷ 
درصد مربوط به واحدکاری قاسملو و کمترین با مقدار ۲۷/۱۷ درصد مربوط به واحد کاری بابارود می‌باشد 
«rani, Abghari & Rasouli, 2024)‏ شکل )1( موقعیت حوزه آبخیز باراندوزچای را در سطح کشور و 


استان آذربایجان غربی نمایش می‌دهد. 


دریاچه ارومیه mE‏ 


Urmia Lake 
| حوضه آبریز دریاچه ارومیه|‎ 
Urmia Lake Watershed 
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Barandozchay Watershed 


شکل ۱- موقعیت جغرافیایی حوزه آبخیز باراندوزچای- دریاچه ارومیه 


Fig.1. Geographical location of the Barandozchay Watershed - Urmia Lake 


۱۶۰۳ جغرافیا و مخاطرات محیطی. جلد ۰۱۳ شماره ۳ پائیز‎ 4 pi VF 


5351-33 م 

در پژوهش حاضر از تصاویر ماهواره‌ای L Sentinel2‏ تفکیک مکانی E‏ متری و دارا بودن محد‌وده طبف‌های 
متنوع (۱۳ باند برای (Sentinel2‏ استفاده گردید. اطلاعات تصاویر مورداستفاده در این مطالعه در حدول (۱) 
و (Y)‏ ارائه شده است. 


جدول ۱- اطلاعات ماهواره Sentinel2‏ 


Table 1- Sentinel 2 satellite information 


باند وضوح طول‌موج مرکزی شرح 
Description Central wavelength Resolution band‏ 
آبی فوق‌العاده (سا sí‏ 
nm 60m BI‏ 443 بی فوق‌العاده (ساحلی و اثروسل) 


Ultra Blue (Coastal and Aerosol) 


[ 
= 490 nm 10m B2 
Blue 
ee 560 nm 10m B3 
Green 
فک‎ 665 nm 10m B4 
Red 
(VISIR) و تزدیک‎ uf e مادون‌قرمز‎ 705 nm 20m B5 
Visible and Near Infrared (VNIR) 
(VNIR) مادون‌قرمز مرئی و نزدیک‎ 740 nm 20m B6 
Visible and Near Infrared (VNIR) 
C(VNIR) و نزدیک‎ als poste 783 nm 20m B7 
Visible and Near Infrared (VNIR) 
(VNIR) مادون‌قرمز مرئی و نزدیک‎ 842 nm 10m B8 
Visible and Near Infrared (VNIR) 
(VNIR) مادون‌قرمز مرثی و تزدیک‎ 865 nm 20m B8a 
Visible and Near Infrared (VNIR) 
SWE) ی‎ 940 nm 60m B9 
Short Wave Infrared (SWIR) 
Short Wave Infrared (SWIR) 
(SWIR) موج‌کوتاهمادون‌ترمز‎ 1375 nm 60m B10 
۱۳۱۳۲۹ pple مۇج كرتا‎ 1610 nm 20m B11 
Short Wave Infrared (SWIR) 
(SWIR) موج‌کوتاه مادون‌قرمز‎ 2190 nm 20m BI2 


Short Wave Infrared (SWIR) 


کوهانی و بهمنش. تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... ۳۷۵ 


حدول ES‏ اطلاعات FE‏ ماهواره‌ای مورداستفاده در مطالعه تغییر و پیش‌بینی کاربری 


Table 2- Information of satellite images used in the study of change and prediction of land use 


نام سنجنده تاریخ کسب داده قدرت تفکیک مکانی 
Spatial resolution Date of data acquisition The name of the gauge‏ 
سنتینل Y‏ 6 معادل 2016.04.26 چندطیفی 10 20 و 60 متر 
Multispectral 10, 20 and 60 meters 2016.04.26 Equivalent 2016.07.16 Sentinel2‏ 
سنتینل Y‏ 0 معادل 1401.03.30 چندطیفی 10 20 و 60 متر 
Multispectral 10, 20 and 60 meters 1401.03.30 Equivalent 2022.06.20 Sentinel2‏ 


۲-۲-۸ - پیش پردازش FE‏ ماهواره‌ای: در al> „o‏ پیش‌پردازش تبدیل فرمت تصاویر» برش "E‏ و 


موزائیک‌سازی در محیط نرم‌افزارهای (ERDAS Imagine and ENVI)‏ انجام گردید. 


۲-۲-۲- در مرحله پردازش: ورود داده‌ها و تصاویر. ایجاد پروژه‌های مورد نیاز» سگمنت‌سازی چند متغیره. 
ویراستاری شی‌های تصویری, طبقه‌بندی اراضی زمین» تولید شاخص‌های طیفی, برآورد دقت و صحت نتایج» 
فیلتر گذاری و ویراستاری ارسال لایه‌های رقومی و طبقه‌بندی نظارت‌شده به روش‌های NNC‏ تولید نقشه‌های 
پوشش زمین و کاربری اراضی در le‏ نرم‌افزار 9.5 eCognition Developer‏ انجام گردید. در ادامه بعد از 
فراخوانی تصوير و ایجاد ترکیبی باندی مناسب. لازم به Sd‏ است که در پژوهش حاضر از باندهای YEA‏ 
و ۱۱ استفاده گردید. بعد از ایجاد ترکیب GUL‏ سگمنت‌سازی (قطعه‌بندی) انجام شد. سگمنت‌سازی به 
مفهوم جداسازی پیکسل‌های تصویری در باندهای مختلف تصاویر ماهواره‌ای و یک پارچه‌سازی مجدد آن بر 
مبنای پارامترهای مختلف نظیر شباهت‌های بازتابشی. هندسی, تراکمی» فشردگی, ارتباط ذاتی آن‌ها باهمدیگر 
از طریق اعمال الگوریتم‌های تقسیم‌بندی متعدد تعریف می‌گردد )2020 (Nazar Neghad et al,‏ نرم‌افزار 
eCognition Developer 9.5‏ جهت طبقه‌بندی راه‌های مختلفی در پیش روی کاربر قرار می‌دهد. 
طبقه‌بندی‌هایی که غالبا در این زمینه و با استفاده از این نرم‌افزار انجام گرفته به‌صورت مراحل سگمنت‌سازی» 
نمونه‌برداری و طبقه‌بندی است. ازآنجایی که در alm‏ طبقه‌بندی شیءگرا تعریف شی‌ها یا سگمنت‌ها است» 


به همین منظور برای سکمنت‌سازی JN «0 E‏ و ۰/۷ به ترتیب برای پارامترهای مقیاس. شکل 3 


i ۳۷۶‏ 4 جغرافیا و مخاطرات محیطی. جلد ۰۱۳ شماره ۳ پائیز ۱۶۰۳ 


فشردگی مطلوب بودن در نظر گرفته شد. جهت تفکیک up‏ پوشش گیاهی از شاخص (NDWI!)‏ و 
(NDVI?)‏ استفاده شد. رابطه شاخحص‌ها برای ماهواره Sentinel2‏ در روابط ۱ و Y‏ نشان داده شده است 
.(Moharrami & Derafshi, 2015)‏ 


NDVI Sentinel2 = (Near Infraredb8 — REDb4) /(Near Infraredb8 + REDb4) Q) 


NDWI Sentinel2 = (Greenb3 — Near Infraredb8) /(Greenb3 + Near Infraredb8) (Y) 


پس از ایجاد شاخحص‌هاء با توجه به هدف پژوهش که آشکارسازی تغییرات کاربری اراضی است. کلاس‌های 
کاربری در نظر گرفته‌شده شامل کاربری پهنه‌آبی. مسکونی. مرتع. نمک حاشیه دریاچه ارومیه باغات و 
اراضی cel‏ اراضی دیم و خاک است. در طبقه‌بندی تصوير با روش شیءگرا از الگوریتم نزدیک‌ترین همسایه 
NNC‏ استفاده شد و متناسب با Ghee per‏ بصری و به کمک شاخص‌ها نمونه برای هرکدام از کلاس‌ها 
انتخاب و شرایط طبقه‌بندی برای هر یک از کلاس‌های کاربری اراضی تعریف شد. نقشه‌های کاربری اراضی 
برای Y‏ دوره با روش شیءگرا استخراج dos S‏ و مساحت کاربری‌ها طبقه‌بندی شده محاسبه گردید. 

۳-۲-۲- محاسبه سنجه‌های سیمای سرزمین 

یکی از ابزارهای محاسبه سنجه‌های سیمای سرزمین» نرم‌افزار Fragstats‏ می‌باشد که استفاده از آن به d^‏ 
قابلیت بیش‌تر و محاسبه تعداد بیش‌تری از سنجه‌های سیمای سرزمین معمول‌تر است )2002 (McGarigal,‏ 
نرم‌افزار 125 مجموعه کاملی از سنجه‌های سیمای سرزمین می‌باشد. این برنامه محدودیتی در مقیاس 
(پهنه و بزرگ‌نمایی) نداشته و برای تجزیه‌وتحلیل الگوی مکانی و سنجه‌های مختلف لکه‌های تشکیل‌دهنده 
سیمای سرزمین در محیط‌های ناهمگن و شرایط متفاوت مناسب است. به‌ازای هر ورودی در نرم‌افزار 
a. Fragstats‏ فایل خروجی Class Patch‏ و Landscape‏ در محیط Excel‏ ذخیره شده و قابل مشاهده و 
تجزیه‌وتحلیل است. آنچه که سنجه‌های سیمای سرزمین را کاربردی‌تر ساخته است. امکان مقایسه سیمای 


سرزمین‌های مختلف با استفاده از یک روش مشابه» ارزیابی در یک سیمای سرزمین yee‏ برای دفعات متعدد 


! Normalized Difference Water Index 


? Normalized Difference Vegetation Index 


کوهانی و بهمنش. تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... ۳۷۷ 


و مقایسه سیمای سرزمین‌های مشخص با گزینه‌های متفاوت است )2012 (McGarigal, Cushman, & Ene,‏ در 


شکل (۲) فرآیند IS‏ انجام پژوهش ارائه شده است. 


دریافت و مدیریت تصاویر ماهواره‌ای 56011061-2 


Acquisition and management of satellite 
imagery Sentinel-2 


Jil‏ تصاویر به محیط نرم افزار 


eCognition Developer 9.5 Pre-processing of images 
Transfer of images to the software (ERDAS Imagine & ENVI) 
eCognition Developer 9.5 


پیش - پردازش تصاویر 


ل 


تولید Lb cls rls‏ اعمال شروط آستانه‌گذاری و فیلترهای قانونمند 
qas Mie‏ قطعه‌بندی چند مقیاسی- 
Production of spectral‏ تطبیقی و شئ گرا 
A : indices‏ ; تفکیکی تصاویر ماهواره‌ای 
Í . Applying conditions, thresholding, and NDVI and NDWI‏ 
Multi-scale segmentation rule-based filters adaptive and object-‏ 
Satellite imagery oriented‏ 


طبقه‌بندی نظارت‌شده به روش NNC‏ 


سنجه‌های سیمای سرزمین 


NNC Supervised classification method 


نرم افزا 
"WT‏ |[ نرم افزار 
0 ` تولید نقشه‌های کاربری اراضی Landscape indicators‏ 
Software Fragstats Change Detection‏ 


ArcGIS 


تعداد (NP) aS‏ شاحص بزرگ‌ترین لکه (LPD‏ شاحص شکل لکه (LSD‏ تراکم حاشیه ED)‏ حاشیه کل 
(TE)‏ تراکم لکه (PD)‏ میانگین اندازه لکه (AREA -MN)‏ میانگین فاصله نزدیک‌ترین همسایه اقلیدسی 
(ENN-MN)‏ شاخص گسستگی سيما (DIVISION)‏ میانگین شاخص شکل «SHAPE -MN) «SJ‏ 
شاخحص تکه‌شدگی (SPLIT)‏ تراکم غنای galė (PRD) aS‏ تنوع شانون (SHDID)‏ شاخص اصلاح 
شده تنوع سیمپسون (MSIDD‏ شاخص پیوستگی لکه (COHESION)‏ شاخص تجمع (AT)‏ اندازه شبکه 
تاثیر گذار (MESH)‏ شاخص پراکندگی (III)‏ درصد مجاورت‌های مشابه (PLADJ)‏ شاخحص سرایت 
(CONTAG)‏ میانگین au‏ همسایگی (CONTIG MN)‏ میانگین شاخص ابعاد فراکتال 
(FRAC_MN)‏ 


تجزیه و تحلیل ساختار چشم‌انداز در دو سطح کلاس و سیما 
Analysis of landscape structure at two levels: class and image‏ 


شکل ۲- نمودار جریانی مراحل انجام پژوهش 


Fig. 2. Flow chart of research steps 


۱۶۰۳ شماره ۲ پائیز‎ OY جغرافیا و مخاطرات محیطی» جلد‎ 4 pia YVA 


۳- نتایج و بحث 

پردازش تصاویر به روش‌های شیءگرا پیکسل‌های تصویر ماهواره‌ای (تشکیل‌دهنده باندهای تصویری) با 
اعمال الگوریتم‌های (Segmentation) T TURO Ld‏ و با توجه به میزان بازتابش» شکل. تراکم و ple‏ 
ارزش‌های زمین-آماری به شی‌های تصویری تبدیل گردید با اعمال روش‌های ایندکس‌سازی: طبقه‌بندی 
تصاویر انجام گردید. در این پژوهش از ترکیب A . ۳ ۲ cloth‏ و ۱۱ تصویر Sentinel2‏ استفاده شده 


است. نتایج ترکیب باندی تصویر» سگمنت‌سازی و طبقه‌بندی تصویر در شکل (۳) ارائه شده است. 


Image band composition 
تصویر‎ sth ترکیب‎ 


Segmentation 
سگمنت سازی (قطعه بندی تصویر)‎ 


Image classification 
طبقه‌بندی تصویر‎ 


شکل ۳- ترکیب باندی. سگمنت سازی و طبقه‌بندی تصویر 


Fig. 3. Band composition, segmentation and image classification 


بعد از طبقه‌بندی تصاویر به روش NNC‏ در هفت کاربری MIT‏ خاک باغات و اراضی "PI‏ اراضی دیم» 
شتا حت ھر es‏ ماه کد کل GY‏ ها eas‏ اراک ال ار PLS VOTE‏ مک 


طبق نتایج برای سال ۲۰۱ مساحت پهنه‌آبی» خاک. باغات و اراضی آبی. اراضی se‏ نمک حاشیه دریاچه 


کوهانی و بهمنش» تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... Yva‏ 


ارومبه. ^( و مسکونی به ترتیب ۵ ٩ E ۵۵۰/۷۲ ۰۰/۰۱ (40/0Y ۰۲۲۱/۱۶ ATIVE‏ کبلومترمربع به 
دست آمد. همچنین در سال ۲۰۲۲ تغییرات مساحت برای کاربری‌های پهنه آب» خاک باغات و اراضی cel‏ 
اراضی qe?‏ نمک حاشیه دریاحه Er cane gl‏ و مسکونی به —- ۰ ۳۳۸/۷۶ ۰۶۳۷/۶۸ ۰۱۱۶/۱۸ ۱/۰۱۵ 


۲۱ و 10/0۰ کیلومترمربع بود. 


44°50'0"E 45°0'0"E 45?10'0"E 44°50'0"E 45°0'0"E 45°10'0"E 45?20'0"1 


37°24'0"N 
37°24'0"N 


37°13'30"N 
37°13'30"N 


44°50'0"E 45°0'0"E 45?10'0"E 44°50'0"E 45°0'0"E 45°10'0"E 


دیم © Rainfed‏ 
مناطق مكينى identialareas‏ 
مرتع ۱0 Rangeland‏ 
نمک حاشیه دریاچه ارومیه © lake Urmia‏ 91 وی 
خاک لا Sol‏ 
یاغات و اراضی d irrigated lands B.‏ 


Water Bodies # jÎ Ste 


شکل ۶- نقشه کاربری اراضی برای سال‌های موردبررسی در حوزه آبخیز باراندوزچای-دریاچه ارومیه 
Fig. 4. Land use map for the studied years in the Baranduzchay catchment area‏ 
بعد از تهیه نقشه‌های کاربری اراضی. آشکارسازی تغییرات برای درک تغییرات منطقه» گسترش و کاهش 


مساحت کلاس‌ها در سال‌های ۹ و ۲۰۲۲ به کیلومترمربع و درصد در حدول (Y)‏ ارائه Cou ados‏ 


۱۶۰۳ شماره ۲ پائیز‎ OY جغرافیا و مخاطرات محیطی» جلد‎ 4 pia YA- 


جدول ۳- مساحت کاربری اراضی به کیلومترمربع و درصد 


Table 3- land use area in square kilometers and percentage 


سال 
2016 2022 
Year‏ 
E‏ کیلومترمریع کیلومترمریع TM‏ 
square square‏ 
Percent kilometer Percent Kilometer Land Use‏ 
ديم 
95.51 7.02 114.18 8.39 
Rainfed‏ 
pm‏ 
10.06 0.74 25.50 1.87 
Residential areas‏ 
eo‏ 
550.73 40.45 412.91 30.33 
pasture‏ 
نمک حاشیه دریاچه ارومیه 
0.01 0 0.05 0 
Salt on the edge of Lake‏ 
Urmia‏ 
خاک 
482.74 35.46 338.74 24.88 
Soil‏ 
SEL‏ و اراضی آبی 
221.14 16.24 467.48 34.33 
Gardens and water lands‏ 
پهنه‌ی آب 
1.35 0.1 2.7 0.2 


water area 


کوهانی و بهمنش, تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... YA)‏ 


ORamfed نیم‎ 

۷ residential areas مناطق مسکوتی‎ 
B pasture مرتع‎ 

B Lakeside salt تمکزار‎ 

7 Soils 


8i Gardens and water lands si یاغات و اراضی‎ 


۷ water area پینه آبی‎ 


شکل ۵- روند تغییرات کاربری اراضی به روش شی گرا 
Fig. 5. The trend of land use changes by object-oriented method‏ 


بر اساس نتایج حاصل از روند تغییرات شکل O‏ در این پژوهش» بیش‌ترین تغییرات کاربری مربوط به 
تبدیل کاربری مرتع و SE‏ به کل اراضی آبی و EL‏ دیم و مسکونی را نشان می‌دهد که با نتایج 
Birhanu, Masih, van der Zaag, Nyssen & Cai, 2019))‏ در پژوهش بر )$9( حوضه PR‏ گومارا' اتیوپی 
نشان دادند در مدت YA‏ سال در اتیوپی مساحت مراتع کاهش aL‏ و در مقابل مساحت اراضی کشت شده 


افزایش یافته و همین امر باعث کاهش حجم آب‌های جاری حوضه آبریز شده است مطابقت دارد. 


تجزیه و تحلیل سنجه‌ها در سطح کلاس 

ae‏ حاصل از ارزیابی تغییرات مکانی سنجه‌های سیمای سرزمین حوزه آبخیز باراندوزچای در CU JRE‏ نشان 
oala‏ شده است. برای اندازه‌گیری میزان گسستگی سیمای سرزمین تعداد لکه (NP)‏ به کار می‌رود و تعداد 
لکه‌های یک کلاس معین را OUS‏ می‌دهد. مقدار آن بزرگ‌تر یا مساوی یک است. وقتی سیمای سرزمین فقط از 
یک fas asd‏ شده th‏ این la‏ پرابر وال است )2015 Riya & Peghhi,‏ برای سرض مور دمطالع ه 


بیش‌ترین تعداد لکه در سال ۱۳۹۵ مربوط به مرتع (۲۷۳۱) و کم‌ترین آن برای پهنه آبی OV‏ و برای سال ۱۶۰۱ 


1 Gomara 


۱۶۰۳ شماره ۲ پائیز‎ Y جغرافیا و مخاطرات محیطی» جلد‎ 4 pia YAY 


پیش‌ترین تعداد لکه (NP)‏ مربوط به مرتع 070( و کم‌ترین آن برای پهنه آبی CW)‏ به دست آمد. 

سنجه تراکم لکه (PD)‏ برای اندازه‌گیری پیوستگی کاربری‌ها و میزان اتصال کاربرد دارد. پیوستگی بین لکه‌ها 
(Cohesion)‏ به‌صورت معکوس با سنجه تراکم لکه ارتباط دارد. به‌این‌ترتیب با افزايش تراکم AS‏ پیوستگی 
سیمای سرزمین کاهش و لکه‌ها منظم‌تر و کوچک‌تر می‌شوند(2020 (Nazar Neghad et al.,‏ بیش ترین تراکم 
450 در سال ۱۳۹۵ مربوط به مرتع (IY)‏ و کم‌ترین آن مربوط به پهنه‌آبی (۰/۰۲) است. برای سال ۱۶۰۱ 
بیش‌ترین تراکم لکه مربوط به مرتع CINE)‏ و کم‌ترین OF‏ مربوط به پهنه‌آبی (۰/۰۵) به دست آمد. نتایج 
پژوهش در ارزیابی الگوی کاربری اراضی و تغییرات پیوستگی نشان‌دهنده افزایش تراکم و تعداد لکه‌ها بود. 
این نتیجه با پافته‌های نظرنژاد و همکاران )2020 (Nazar Neghad et al.,‏ مبنی بر این که افزایش تراکم 3 
تعداد لکه‌هاء افزایش ناهمگنی در کاربری اراضی را نشان می‌دهد. مطابقت دارد. حداکثر مقدار شاخص 
بزرگ‌ترین لکه (LPD‏ در سال ۱۳۹۵ مربوط به (Dg‏ و حداقل مقدار آن مربوط به پهنه آبی (۰/۰۱) 
است و در سال ۱۶۰۱ حداکثر مقدار شاحص بزرگ‌ترین لکه مربوط به کاربری CEL‏ و اراضی آبی )7( و 
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شکل — مقادیر سنجه‌های سیمای سرزمین در سطح کلاس در سال‌های موردبررسی 
Fig. 6. Values of land surface parameters at the class level in the years under review‏ 
در سطح سیمای سرزمین حاشیه کل (TE)‏ اندازه مطلق مجموع حاشیه از طول حاشیه کل یک نوع لکه خاص 
است )2012 (MeGarigal at aL,‏ با افزایش حاشیه کل» پیوستگی و اتصال سیمای سرزمین کاهش و لکه‌ها 
SS‏ و منظم‌تر می‌شوند. اگر هیچ حاشیه‌ای در سیمای سرزمین وجود نداشته باشد مقدار عددی این شاخص 
زمانی برابر صفر خواهد بود. پعنی زمانی که کل سیمای سرزمین و مرز ol‏ دارای یک لکه باشد McGarigal,‏ 
افزایش تراکم (ED) atl‏ تکه‌شدگی در آن زیرحوضه افزایش می‌یابد. مقدار عددی این شاخص برابر صفر 
خواهد بود. زمانی که هیچ حاشیه‌ای در سیمای سرزمین وجود نداشته باشد. یعنی زمانی‌که کل سیمای سرزمین و 


مرز سیمای سرزمین. شامل یک لکه باشد )2002 (M cGarigal,‏ بیش‌ترین تراکم حاشیه (ED)‏ در سال ۱۳۹۵ 


مربوط به مرتع EN‏ و کم‌ترین مقدار OF‏ مربوط به کاربری پهنه آبی (۰/۱۳) است. هم‌چنین» بیش‌ترین مقدار 


کوهانی و qup‏ تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... YAO‏ 


تراکم حاشیه (ED)‏ در سال ۱۶۰۱ مربوط به کاربری مرتع CUD)‏ و حداقل مقدار آن مربوط به پهنه‌آبی برابر با 
(۰/۳۸) به dal duds‏ 

حداکثر مقدار شاخص شکل لکه (LSD)‏ در سال ۱۳۹۵ مربوط به کاربری مرتع (۱/۱۲) و حداقل مقدار آن 
مربوط به کاربری مسکونی (۱/۰۱) و حداکثر مقدار شاحص شکل لکه در سال ۱۶۰۱ مربوط به کاربری‌های 
باغات و اراضی آبی و مرتع (۱/۱) و حداقل مقدار شاحص شکل لکه به ترتیب مربوط به کاربری‌های مسکونی 
و پهنه ul‏ (۱/۰۱) است. 

میانگین اندازه لکه (AREA-MN)‏ یک نوع ماتریس سیمای سرزمین است که بر اساس ویژگی میانگین لکه 
ارائه شده است )2002 ((McGarigal,‏ حداکثر مقدار میانگین اندازه لکه در هرسال ۱۳۹۵ مربوط به کاربری 
خاک  (Y0/££)‏ حداقل مقدار آن مربوط به پهنه‌آبی (4/۲۶) است و در سال ۱۶۰۱ حداکثر مقدار میانگین 
اندازه لکه مربوط به کاربری باغات و اراضی آبی (۲۵/۵۹) و حداقل مقدار آن مربوط به پهنه آبی (۳/۸۲) 
است. معیار نزدیک‌ترین فاصله همسایگی (ENN, MN)‏ با استفاده از هندسه ساده اقلیدسی به‌عنوان 
کوتاه‌ترین فاصله مستقیم بین نزدیک‌ترین همسایه و لکه مرکزی آن کلاس تعریف شده است ) (McGarigal‏ 
at al., 2‏ این معیار با افزایش بی‌نظمی شکل, بزرگ‌تر می‌شود و برای یک لکه مربع شکل برابر با یک 
است 2015 ((Kiyani & Feghhi,‏ حداکثر مقدار میانگین نزدیک‌ترین فاصله همسایگی در سال ۱۳۹۵ مربوط 
به پهنه‌آبی L)‏ مقدار عددی ۱۷۹۵/۹) و حداقل مقدار آن مربوط به کاربری مرتع (با مقدار عددی ۲۸۷/۱۸) 
است. هم‌چنین حداکثر مقدار ميانگین نزدیک‌ترین فاصله همسایگی در سال‌های ۱۶۰۱ مربوط به کاربری پهنه 
آپی EAST‏ و حداقل آن مربوط به کاربری مرتع (۳۰۲/۵۸) است. حداکثر مقدار شاخص میانگین شکل 
(SHAPE-MN) «SJ‏ برای حوزه آبخیز باران‌دوزچای در سال ۱۳۹۵ مربوط به کاربری مرتع (۱/۱۲) و 
حداقل مقدار آن مربوط به کاربری مسکونی (۱/۰۱) است. برای سال ۱۶۰۱ نیز حداکثر مقدار شاخص 
میانگین شکل لکه مربوط به کاربری باغات و اراضی آبی و مرتع (۱/۱) و حداقل مقدار آن مربوط به کاربری 
پهنه آبی و مسکونی (۱/۰۱) است. معیارهای تکه‌شدگی (SPLIT)‏ و گسستگی سیما (DIVISION)‏ به 
توصیف دخالت انسان‌ها و میزان انسجام اشاره دارند )2002 (McGarigal,‏ تکه‌شدگی (SPLIT)‏ بر اساس 


توزیع لکه‌های تجمعی است و به‌عنوان تعداد لکه‌هایی با اندازه ثابت زمانی یا شبکه موثر تفسیر می‌شود 


۱۶۰۳ جغرافیا و مخاطرات محیطی. جلد ۰۱۳ شماره ۳ پائیز‎ 4 pi TAF 


(McGarigal at al., (2012‏ این معیار حداکثر مقدار را دارد؛ وقتی که سیمای سرزمین به لکه‌های کوک 
تقسیم می‌شود.یعنی زمانی که هر سلول یک لکه جداگانه باشد )2002 (McGarigal,‏ حداکثر مقدار تکه‌شدگی 
در سال ۱۳۹۵ مربوط به پهنه‌آبی (۲۸۶۷۰۰۵/۲) و حداقل OF‏ مربوط به کاربری مرتع (WT)‏ است. 
هم‌چنین» برای سال ۱۶۰۱ حداکثر مقدار تکه‌شدگی مربوط به پهنه‌آبی (۱۶۳۵۵۱۲۵/۲) و حداقل Ul‏ مربوط 
به کاربری باغات و اراضی - (EY)‏ به دست آمد )2014 (Su, Naqui, Mata-Contreras & Martín,‏ در 
بررسی تکه‌شدگی (SPLIT)‏ سیمای سرزمین از سنجه‌های حاشیه کل (TE)‏ تعداد لکه (NP)‏ شکل سیمای 
سرزمین» میانگین سطح لکه (AREA -MN)‏ بزرگ‌ترین لکه (LPD‏ میانگین شاخص شکل S)‏ 
(SHAPE -MN)‏ استفاده شد. نتایج نشان داد که زراعت می‌تواند انجام شود. مقدار گسستگی سیما 
(DIVISION)‏ در حوزه آبخیز باراندوزچای در سال ۱۳۹۵ بین ۱-۰/۹۹٩۹‏ و در سال ۱۶۰۱ بین ۱-۰/۹۸ متغیر 
بوده است. این نتایج با نظرنژاد و همکاران )2020 (Nazar Neghad et al.,‏ که در سال‌های 1٩۹۰٩ AAA‏ 
۷۱ و ۲۰۰۲ به این نتیجه رسیدند تنوع در سیمای سرزمین کاهش OL‏ و نواحی علفزار شهری. مرتع و 
شالیزار در داخحل لکه‌های بزرگ متمرکز شده‌اند و جداشدگی مربوط به پهنه آبی بوده است. مطابقت دارد. 
حداکثر مقدار شاخص تجمع (AD‏ در سال ۱۳۹۵ مربوط به کاربری باغات و اراضی آبی )۵۹/٩۱(‏ و حداقل 
مقدار آن مربوط به کاربری مسکونی (۱۳/۹۵) است. در سال ۱۶۰۱ نیز حداکثر مقدار شاخص تجمع مربوط 
به CEL‏ و اراضی آبی (11/1۳) و حداقل مقدار ol‏ مربوط به پهنه‌آبی MOY)‏ به دست آمد. شاخص 
پیوستگی (COHESION) aS‏ ارتباط فیزیکی لکه‌های مربوط به یک کاربری را بیان می‌کند. مقدار این 
شاخص با افزایش پیوستگی لکه‌ها نیز افزایش خواهد یافت. ) Alaei, Mustafazadeh, Esmaaliouri, Sharari‏ 
Hazbavi, 0‏ &( حداکثر مقدار شاخص پیوستگی لکه در سال ۱۳۹۵ مربوط به کاربری SLE‏ )46( و 
حداقل مقدار آن مربوط به پهنه‌آبی (۱۷/۳) است. در سال ۱۶۰۱ حداکثر مقدار این شاخص مربوط به کاربری 
باغات و اراضی آبی (4۵/۲) و حداقل مقدار OF‏ مربوط به پهنه‌آبی QVO‏ به دست آمد. 

تجزیه و تحلیل سنجه‌ها در سطح سیما 

نتایج حاصل از ۲۲ سنجه حوزه آبخیز باراندوزچای جدول (۳) oala OLS‏ شده است. نتایج مشخصه‌های 


آماری این سنجه‌ها می‌توان نتیجه‌گیری کرد که حداکثر تراکم لکه مربوط به سال ۱۶۰۱ و حداقل مقدار 


کوهانی و qup‏ تحلیل ارتباط بین تغییرات کاربری‌اراضی و سنجه‌های سیمای سرزمین ... YAV‏ 


مربوط به سال ۱۳۹۵ بوده است. وجود تراکم بالای لکه در حوزه آبخیز نشان‌دهنده تخریب کاربری اراضی 
است که باعث افزایش تکه‌شدگی نیز شده است. پایین بودن سنجه میانگین اندازه ASS‏ و بزرگ‌ترین اندازه 
لکه به معنی تکه‌شدگی بیش‌تر سیمای سرزمین است. مقدار شاخص بزرگ‌ترین لکه در سال ۱۳۹۵ با مقدار 
عددی ٩/٩‏ و کم‌ترین OF‏ در سال ۱۶۰۱ با مقدار ۹/۱ به دست آمد. با افزایش و کاهش مقدار عددی شاخحص 
بزرگ‌ترین لکه به ترتیب میزان تکه‌شدگی سیمای سرزمین کاهش و افزايش می‌یابد. ازاین‌ری در سال ۱۳۹۵ 
Ol jes‏ تکه‌شدگی کم‌ترین مقدار و در سال ۱۶۰۱ بیش‌ترین مقدار می‌باشد. هر چه تراکم حاشیه (ED)‏ بیش تر 
باشد اتصال و پیوستگی سیمای سرزمین بیش‌تر می‌شود؛ بنابراین با توجه به جدول Y)‏ حداکثر اتصال بین 
لکه‌ها در سال ۱۶۰۱ بوده و در سال ۱۳۹۵ حداقل تراکم حاشیه (ED)‏ و درنتیجه بیش‌ترین تکه‌شدگی 
مشاهده می‌شود. سنجه شکل لکه (LSD)‏ نسبت محبط لکه به حداقل محیط ممکن برای یک لکه با حداکثر 
تجمع را OLE‏ می‌دهد. زمانی که مقدار عدد آن برابر یک باشد. یعنی سیمای سرزمین شامل یک لکه با 
حداکثر فشردگی و تقریباً مربع شکل است» وقتی که لکه پراکنده‌تره شکل آن پیچیده‌تر و مرز بی‌نظم‌تر شود 
مقادیر آن بدون محدودیت افزایش می‌یابد. با توجه به جدول (۳)» چون در اکثر زیر حوضه‌های آبریز بیش تر 
از یک می‌باشد می‌توان نتیجه گرفت شکل لکه‌ها در تمامی سال‌های موردبررسی نامنظم می‌باشد. میانگین 
شاحص JS‏ لکه (SHAPE MN)‏ از نسبت محیط لکه به بیش‌ترین فشردگی یا کم‌ترین محیط ممکن 
برای لکه به دست می‌آید. این سنجه برای لکه مربع شکل برابر یک است و با نامنظم شدن شکل, بیش‌تر 
می‌شود. با توجه به جدول O)‏ روند تغییرات SHAPE MN‏ برای سال‌های موردبررسی روند کاهشی کمی 
داشته است» به‌طور ی که از مقدار عددی ۱/۰۹ در سال ۱۳۹۵ به مقدار ۱/۰۸ در سال ۱۶۰۱ رسیده است. 
حداکثر کل حاشیه (TE)‏ 14۰1۸۰۰ برای سال ۱۳۹۵ و برای سال ۱۶۰۱ مقدار ۷۳۹۹۲7۰ ارزیابی شد. طبق 
تغییرات این سنجه روند به‌صورت افزایشی بوده است. با افزايش حاشیه کل پیوستگی و اتصال سیمای 
سرزمین کاهش و لکه‌ها منظم‌تر و کوچک می‌شوند )2012 (McGarigal at al.,‏ شاخص پیوستگی لکه‌ها 
(COHESION)‏ ارتباط فیزیکی لکه‌های مربوط به کاربری‌ها را OUS‏ می‌دهد. با افزایش پیوستکی لکه‌ها 
مقدار این شاحص نیز افزايش خواهد یافت. با توجه به نتایجم جدول polis (E)‏ به دست آمده برای تمام 


بخش‌های حوزه آبخیز کاهشی است. بنابراین می‌توان نتبجه گرفت که پیوستگی laa S‏ در تمامی سال‌های 
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موردبررسی روند کاهشی دارد. سنجه‌های شانون و سیمپسیون از مهم‌ترین سنجه‌ها برای مطالعه تنوع اجزای 
تشکیل‌دهنده سیمای سرزمین هستند. شاخص شانون. تنوع را در هر لکه به‌طور نسبی محاسبه می‌کند. این 
سنجه‌ها زمانی برابر با صفر هستند که در سیمای سرزمین فقط یک AS)‏ وجود داشته باشد. و زمانی مقدار 
آن‌ها یک خواهد بود که تعداد لکه افزایش می‌یابد و توزیع منطقه متناسب با انواع لکه بیش‌تر شود 
(McGarigal at al., 2012)‏ سنجه شانون (SHDI)‏ نسبت به نادرترین لکه‌ها حساس است. زمانی میزان تنوع 
افزایش می‌یابد که میزان سنجه‌های تنوع به یک نزدیک شود. هرچه این مقدار به صفر نزدیک شود 
از ميزان تنوع کاسته می‌شود. 
جدول -٤‏ مقادیر سنجه‌های سیمای سرزمین در سطح سیما در سال‌های موردبررسی 


Table 4- The values of land surface measurements at the surface level in the years 


under review 


روند تغییرات 
نام فارسی سنجه احتصار 
Trend of 2022 2016‏ 
changes Abbreviated The Persian name of Sanjeh‏ 
تعداد لکه افزایشی 
NP‏ 7162 9858 
increasing Number of Patches‏ 
تراکم SS‏ افزایشی 
PD‏ 5.7 7.24 
increasing Patch Density‏ 
شاخص بزر 4S) n pe‏ کاهشی 
LPI‏ 9.6 9.1 
decreasing largest patch index‏ 
حاشیه کل افزایشی 
TE‏ 6904800 7399260 
increasing Total Edge‏ 
تراکم حاشیه 
افزایشی 
ED Edge Density‏ 50.72 54.35 
increasing‏ 
شاخص شکل لکه gestit‏ 
LSI‏ 48.52 51.88 
increasing Landscape Shape Index‏ 
میانگین اندازه لکه کاهشی 
AREA MN‏ 54.17 13.81 
decreasing Mean patch area‏ 


میانگین شاخص شکل لکه SHAPE MN‏ 1.09 1.08 کاهشی 
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decreasing Mean patch shape index 
cat میانگین شاخص ابعاد فراکتال‎ 
1.01 1.01 FRAC MN 
constant Fractal Dimension Index 


میانگین شاخص همسایگی 


ثابت 
CONTIG MN Contiguity Index‏ 0.16 0.16 
constant‏ 
میانگین فاصله نزدیک‌ترین همسایه 
اقلیدسی کاهشی 
ENN MN‏ 371.53 361.07 
decreasing Mean Euclidean Nearest‏ 
Neighbor Distance‏ 
شا خن رابت کاهشی 
CONTAG‏ 39.61 37.83 
decreasing Contagion index‏ 
درصد مجاورت‌های مشابه کاهشی 
PLADJ xd‏ 53.50 50.23 
decreasing Percentage of similar‏ 
adjacencies‏ 
شاخص پراکندگی افزایشی 
IJI‏ 51.28 56.82 
increasing dispersion index‏ 
شاخص پیوستگی ASS‏ کاهشی 
COHESION‏ 93.36 88.19 
decreasing Patch Cohesion Index‏ 
شاخص گسستگی سیما ثابت 
DIVISION‏ 0.98 0.98 
constant Landscape Division Index‏ 
sls‏ 
اندازه شبکه SE‏ گذار MESH‏ 8 2383.13 
decreasing‏ 
شاخص تکه‌شدگی افزایشی 
SPLIT‏ 42.04 57.13 
increasing Splitting Index‏ 
تراکم غنای لکه ثابت 
PRD‏ 0.01 0.01 
constant Patch Richness Density‏ 
شاخص تنوع شانون 
افزایشی 
SHDI Shannon's Diversity Index‏ 126 137 
increasing‏ 
شاخص اصلاح شده تنوع سیمیسون افزایشی 
ACA RN ne MSIDI Modified Simpson's‏ 
Diversity Index‏ 8 
شاخحص تج کاهشی 
AI S‏ 54 50.71 


decreasing Aggregation index 
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-٤‏ نتیجه گیری 

در پژوهش حاضر به‌منظور بررسی تغییرات کاربری اراضی و سنجه‌های سیمای سرزمین در حوزه آبخیز 
باراندوزچای با لحاظ معیارهای موثر بر آن شامل ۲۲ سنجه در سطح سیما و ۱۳ سنجه در سطح کلاس 
محاسبه شد. سنجه تکه‌شدگی نشان داد که کاربری پهنه‌آبی حداکثر تکه‌شدگی و زراعت دیم حداقل 
تکه‌شدگی را دارا می‌باشند. نتایج OUS‏ داد که با توجه به سنجه‌های تراکم aS)‏ شاخص بزرگ‌ترین aI‏ 
تراکم حاشیه و شاخص شکل لکه کاربری مرتع حداکثر و کاربری زراعت دیم حداقل پیوستگی بین لکه‌ها را 
دارند. در این پژوهش با توجه به نتایج حاصل از سنجه‌های تنوع شانون و سیمسیون در سطح سیما 
مشاهده شد که اگر مقدار این نسبت dos‏ 39 باشد. پراکنش کپه‌ای و اگر Apes erase. la‏ 
درصورتی که کم‌تر از یک باشد به‌صورت یکنواخت توزیع یافته‌اند بنابراین همان‌گونه که مشاهده شد 
این مقدار برای بازه زمانی ۱۳۹۵ و ۱۶۰۱ بیش از یک به دست آمد. پس می‌توان نتیجه گرفت که در حوزه 
seul‏ باراندوزچای پراکنش نکه‌ها به‌صورت یکنواخت می‌باشد. GUS quld‏ داد که افزایش تعداد تکههاه 


افزايش حاشیه و کوچک شدن لکه‌ها باعث تکه‌شدگی و درنتیجه ایجاد لکه‌های جدا از هم در حوزه آبخیز 


باراندوزجای می‌شود. نتایج این پژوهش می تواند 22 — کاربری اراضی. جلوگیری از بهره‌برداری 


بیش ازحد در منطقه موردمطالعه و حفاظت از Of‏ برای نسل‌های آتی مفید باشد. 
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